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CALCULATION OF LAMINAR AND TURBULENT BOUNDARY LAYERS 
FOR TWO-DIMENSIONAL TIME-DEPENDENT FLOWS 

by 

Tuncer Cebeci 
I. SUMMARY 

In recent years several methods have been developed to compute time- 
dependent laminar and turbulent boundary layers. Some of these methods use 
the eddy- viscosity, mixing-length concepts to model the Reynolds stresses 
(for example, see references 1 and 2), and others use the approach advocated 
by Bradshaw (see for example, references 3 and 4). The chief advocates of 
Bradshaw’s method are Nash and Patel. 

A general method for computing laminar and turbulent boundary layers for 

two-dimensional time-dependent flows is presented in this report. The method 

uses an eddy-viscosity formulation to model the Reynolds shear-stress term and 

a very efficient numerical method to solve the governing equations. The model 

f 5) 

developed by Cebeci, has been applied to two-dimensional flows' ' and to three- 
dimensional flows^^*^^ and has been shown to give good results. 

In Section II the governing equations and the eddy-viscosity formulas are 
presented. When physical coordinates are used, the numerical solutions of the 
boundary-1 ayer equations are quite sensitive to the spacings in the x- and 
t-directions. In problems where computation time and storage become important, 
it is desirable to reduce the sensitivity to At- and Ax-spacings. This can 
be done by expressing and by solving the governing equations in appropriate 
transformed coordinates and transformed variables. In Section II we introduce 
such transformations. 

Section III presents a discussion of the numerical method and in Section IV 
the results obtained by the present method for both laminar and turbulent flows 
are discussed. Based on the results given in this report, we find the method 
quite efficient and suitable for solving time-dependent flows for both laminar 
and turbulent boundary layers. 



II. GOVERNING EQUATIONS 


2.1 Boundary-Layer Equations 

The continuity and momentum equations for incompressible, unsteady laminar 
and turbulent boundary layers are: 

Continuity 


3U ^ 9V _ Q 
3X 3y 


( 1 ) 


Momentum 
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( 2 ) 


The boundary conditions for the wall and for the outer edge of the 
boundary layer are: 

y=0 u,v=0 (3a) 

y ” u ^ Ug(x,t) (3b) 

To complete the formulation of the problem, additional boundary conditions 
must be specified on some initial upstream surface, t = T(x) say, normal to 
the x-t-plane. Here we consider the simple case where this surface is made up 


of the two surfaces 
boundary conditions 

t = t^ = const. 

a 

are: 

and X = Xg = const. The required 



and ^ ^ 

u ^ u^(x,y) 

(4a) 

X 

and ^ L 

; u = Ug(t,y) 

(4b) 


Thus at some initial time, conditions are specified everywhere and for subse- 
quent time, conditions are specified on some upstream station. 


2.2 Closure Assumptions for the Reynolds Shear Stress 

The solution of the system given by (1) to (4) requires closure assumptions 
for the Reynolds shear stress - a' v^. In our study we use the eddy viscosity 
concept and define 


—I — r 9U 

-pu V = pe^ ^ 


(5) 


2 



According to the eddy viscosity formulation of ref, 5, the turbulent 
boundary layer is divided into two regions, called inner and outer regions, 
and the eddy viscosity is defined by separate formulas in each region. They 
are: 


m 


(e^)i = j0.4y [1 - exp(-y/A)] 

00 
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3U 

) 

ay 


^^m^i — ^^m^o 

(6a) 

^^m^o — ^^m^i 

(6b) 


0 

In (6a), A is a damping-length constant. For an incompressible flow with no 
mass transfer, it is given by 




+ 1/2 

= u^(l-11.8pp , 


+ 






(7a) 

(7b) 


The parameters a is an "universal" constant equal to 0.0168 for high Reynolds 

number flows, R. > 5000. At low Reynolds number its variation can be approxi- 

0 

mated by a formula given by Cebeci 


a 


a 


1.55 
0 1 + n 


n = 0.55[1 - exp(-0.243z]/^ - 0.298z^)] 


( 8 ) 


where z-j = (Rq/425 — 1). 


Equation (7a) was proposed by Cebeci Por a flqt-plate flow 
(Pj( = 0)» reduces to a form proposed by Van Driest % Cebeci ’s 
extension was made by following Van Driest 's modeling of the viscous sublayer 
to Stokes flow. The shear-wave propagation velocity in Stokes flow was 
assumed to be the friction velocity obtained at a y^- value determined by the 
intersection of the linear law of the wall with the log law of the wall 
(rather than its wall value as was assumed by Van Driest) that is. 


(wv) 


1/2 




1/2 


(9) 
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1 /2 

The friction velocity intersection of the linear law of the 

wall with the log law was determined from the momentum equation approximated 
close to the wall by (no mass transfer) 


dx _ ^ 
dy 3x 


The solution of (10) with y^ = y^/^j u^ = 11.8, 


that is. 


( 10 ) 


1 -n.8p^ (11) 


enables the damping-length constant A to be written in the form given by 
(7a). 


Although this type of approach is highly speculative, and may have little 

theoretical basis (if any), the expression for A obtained in this manner 

(5) 

worked well for turbulent boundary layers with heat and mass transfer^ . 

As an example, we compare the predictions of this approach for an incompres- 
sible flow with no mass transfer. In reference 8, Back, et al , presented some 
results on the generation of turbulence near a smooth boundary in a variable 
free-stream velocity flow by using an analysis based on the extension of 
Einstein and Li's model' . According to the predictions of that model, 
the ratio of the sublayer growth frequency (o to that of a flat-plate, , 

varies with c in the manner shown in Figure 1. The parameter c on which 
the frequency of sublayer growth depends is given by 

5 = 1.21UjP^ (12) 

where ut is the nondimensional velocity at the edge of the sublayer, u /u 

S o T 

and was assumed to be 15.6 in Einstein and Li's analysis as well as in Back's 
analysis. 

According to (9), co is given by 



4 



I 


2.0 

1.6 

1.2 



f.p. 


DECELERATION 
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Figure 1. Comparison of predictions of the present model with those of 
Einstein and Li’s model. 
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We see that with Ug = 15.6, equation (11) can also be written as 

— ^ = 1 - 0.625c (14) 

“f.p. 

The predictions of (14) with those predicted by the extension of Einstein-Li 's 
model are shown in Figure 1. 


Cebeci's extension of Van Driest model to unsteady flows is described in 
reference 10. For completeness, we include it here. For an incompressible 
flow with no mass transfer, integration of (10) between the wall and y^, 
yields 


Noting that. 




T — T * 1 * y 

S W 9X •'s 


9 £ = 

3X 


/ 9U 


e 



(15) 


(16) 
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we can write (15) as 


1 -n.8 (p* + p*) (17) 

i; X 

Here is given by the second expression in (7b) and pt is defined by 

X t 

J. 

e 

Thus (7a) becomes 


1 + + -■'/2 
A = 26vu"' [1 - 11.8(p^ + p^)] 


(19) 


2.3 Transformation of the Governing Equations 

The boundary-layer equations (1) to (2) subject to (3) and (4) can be 
solved when they are expressed either in physical coordinates or in transformed 
coordinates. Each coordinate has its own advantages. In problems where the 
computer storage becomes important, the choice of using transformed coordinates 
becomes necessary, as well as convenient, since the transformed coordinates 
allow large steps to be taken in the x- and t-directions. The reason is 
that the profiles expressed in the transformed coordinates do not change as 
rapidly as they do when they are expressed in physical coordinates. The use 
of transformed coordinates stretches the coordinate normal to the flow and 
takes out much of the variation in boundary- layer thickness for laminar flows. 
In addition, they remove the singularity that the equations in physical 
coordinates have at x = 0 for all time. 


We first use the transformed coordinate n defined by 



and introduce a dimensionless stream function f(x, t, n) defined in 

= >^vxUq f(x,t,n) 


(20) 


( 21 ) 
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Here Uq(x) is some reference velocity and 


u 


ay 



( 22 ) 


With the relations defined by (20) to (22) and with the definition of eddy 
viscosity, (5), it can be shown that the momentum equation can be written as 


Here primes denote differentiation with respect to n and 


f, - U p . X ''“e P - X /„ '’“e . k - i + + 

^ 'IT ’ ■ iraj“ ’ *^3 ■ T' v“e 

0 0 Uq 

The boundary conditions given by (3) become 

n = 0 f = f'=0 


n ->■ n. 


f = u^/u^ 
e 0 


(23) 


(24) 

(25a) 

(25b) 


Similarly the initial conditions given by (4) can be transformed. For some 
problems where the flow starts as laminar, they can be obtained from (23). 
Along the n-t plane where x = 0, (23) reduces to 


ff- _ p(f )^ + p. 


% 


(26) 


Note that, in general, the right-hand side of (26) is not equal to zero. For 
example, although for a flat-plate flow it is equal to zero, for stagnation- 
point flow (u^ = Ax), it is not. For generality we will keep it in that 
form. 


The specification of the initial conditions at t = t^ and x > x, requires 

a — a 

special treatment. Here for simplicity we assume that they are given by the 
steady-state conditions. Then at time t = 0, (23) reduces to 

(bf")‘ + ff" - P(f)2 + Pj = x(f If--f" |i) (27) 


7 


where, for example. 


p - X „ e 

" 77 H)T 


In terms of transformed variables, the eddy viscosity formulas can be 
written as 


(e*)f = 0.16RyV -exp(-y/A)]' 


Here the subscript refers to the boundary-layer edge and 

R ^=-k_y_ n n 4. n+nl/2 

Rx V ’ A 26 ^ + Px)^ 



III. NUMEr?ICAL METHOD 


We use a two-point finite-difference method to: solve the system given by 
(23), (25)-(27). This method, developed by H. B. Keller^^^\ has been applied 


.( 12 ) 


to three-dimensional flows 


to two-dimensional flows by Keller and Cebeci 
by Cebeci Here we apply it to two-dimensional time-dependent flows. 

We first consider (26) and write it in terms of a first-order system of 
partial differential equations. For this purpose we introduce new variables 
u(t,n), v(t,ri) so that (26) can be written as 

f = u (30a) 

u' = V (30b) 

(30c) 


+ P,fv ~ Pu^ + ^ 

1 3 Uq 9t 


Here 


P = P 1 
{^1 2 


We next consider the net rectangle shown in Figure 2. We denote the net 
points by 



Figure 2. Net rectangle for equation (26). 
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t = 0 = t T + k„ n = 1, 2, N (31a) 

0 n n-1 n 

tIq = 0 rij = Tij_i + hj j = 1, 2, J 

Tij = (31b) 

Here the net spacings, k and h. are completely arbitrary, and indeed may 

1 1 J 

have large variations in practical calculations. The variable net spacing h. 

J 

is especially important in turbulent bound ary- layer calculations which are char- 
acterized by large bound ary- layer thicknesses. To get accuracy near the wall, 
small net spacing is required while large spacing can be used away from the 
wall. 


We approximate the quantities (f, u, v) at points (t^, n^) of the 
net by net functions denoted by (f*J, u"!), v!?). We also employ the notation, 

J <J u 

for points and quantities midway between net points and for any net function 



"n-1/2 ■ 2 ^^n ‘‘^n-1^* "j-1/2 

„n-l/2 1 / n . „n-lv „n - j. \ 




(32) 


The difference equations which are to approximate (30) are formulated by 
considering one mesh rectangle as in Figure 2. We approximate (30a, b) using 
centered difference quotients and average about the midpoint 
of the segment P2^4* 


pn _ ^n 

j ^ 




j-1 

" ^j-1/2 

(33a) 


' 3-1 

= v" 

(33b) 


''j-1/2 


Similarly (30c) is approximated by centering about the midpoint 
nj_i /2 of the rectangle PiP2*^3^4* gives: 



py(fv)5.i/2 = T-;^^(33C) 
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where 




^j-1/2 ^“j-V2 '^3' 


The boundary conditions are: 


rv^-L w"-T 
— fr + P-i 

V 

„n-1 / 2vn-l 
-P (uVi/2 


(34) 


f" = 0, u" = 0, u = 


O'T^ 


On 


(35) 


We next consider (27) and again write it in terms of a first-order system 
by introducing u(x,n)» v(x,n). This allows (27) to be written as (30a, b) and 


(bv)' ^p^fv- Pu2 + P3 = x(u|^-v|£) 


(36) 


The net rectangle is the same as before, except now t is replaced by 
X and the net points in the x -direction are represented by 


Xq - 0 , 

The difference equations to (30a,b) are 


^i ^ ^i-1 ’^i* 1=1,2, 


h, “j-1/2 


“3 ~ “3-1 - 
h, ''3-1/2 


Similarly the difference equations for (36) are: 

■ ^ (p] + o,)(fv)1.,/2 - (P’ + “iH“^)]-1/2 


(37) 


(38a) 

(38b) 


(38c) 


n 


Here 


S 


i-1 

j-1/2 






i-1 

j~l/2 





(39a) 


(39b) 


We now consider the general equation (23) and write it in terms of a 
system of first-order equations. We introduce variables u(x,t,n), v(x,t,n) 
so that (23) can be written as 

f = u (40a) 

u’ = V (40b) 

(bv)' .p^fv- Pu2 . P 3 = x(u|^-v|i*l-|^) (40c) 

We now consider the net cube shown in Figure 3. 


n(j) (j,n,i-l) (j,n,i) 



Figure 3. Net cube for the difference equations for (40). 
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For the net points given by (31) and (37), we 'approximate the quantities 
(f, u, v) at points (x^, t^^, nj) of the net by net function denoted by 
(fV”, equations (40a,b) are approximated using centered 

difference quotients and averaged about the midpoint (x^ , t^, '^j-1/2^* 
that is 




''j 




where, for example, 


(bv)j — (b\0^ 1 xi-1/2 




X 


i-1/2 


where, for example. 


■ 1/2 


fUn 

j-1 , ,,i>n 
''j-V2 

(41a) 

u’>'’ 
J-1 ^ 

(41b) 

= i (u1>" + u1>?) 

2 ^ 3 j-r 


1 are to approximate (40c) are; 



-1/2 

-1/2 

“1-1 > - Vi/a^^^i 


, 

"TiipTrK'Vi) 

(42) 

'-1 + .l.(,i.n^^234, 

(43a) 

^ “j-1 ‘ 2 ^“j-1/2 “j-1/2^ 

(43b) 

“j-1 ‘ 2 ^“j-1/2 “j-1/2' 

(43c) 

-1 

(43d) 
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Here by v/e mean 

values of v- at three of the four corners of the face of the box. Intro- 
ducing (43) and similar definitions for other terms into (42), after consid- 
erable algebra, v/e get 

(bv)j - * f'j j “ 

""T [''j-l/2^“j-l/2 * “j-1/2 * “j-1/2 “ “ ''j-l/2*'j-l/2 

“ ^^j-1/2 " ^^i-l'''j-l/2 " ''j-l/2^j-l/2]'' ^Vj-1/2 I " ’’’j-1/2 ^ 


Here 


Bn = 


i-1/2 




T^rTT 


f44) 

(45a) 


j-1/2 ^ "j I 'M'n-1/2^"'j-l/2 ^ '‘^'n-l/2''“ 'j-1/2 

A/n W~1/2 , 1 fll^ ?l7 ^ PT il 

^^^3'n-V2 L j-V2'“j-l/2 ^“i-r ''j-1/2'^0-1/2 

4 26„(u]:]; 5 . 27„.,)j (45b) 


We note that, for simplicity, we have dropped the superscripts 1, n in the 
above equations. 


The boundary conditions given by (25) become 
f„ = 0 , «0 = 0 . “0 = ^ 

The difference equations (33) and (38) for the initial conditions along 
the t-axis, and along the x-axis, respectively, together with the difference 
equations (41), (44), for the general case are imposed for j = 1, 2, J, 
If v/e assume ^ j assume (fj \ 

i 0 ^ j £ J, then (33) and (38) for 1 £ j J and their 

boundary conditions yield an implicit nonlinear algebraic system of 3J + 3 
equations in as many unknov/ns (f*J, u*J, v*J) and (f^, u^-, v^). This system 
can be linearized by using Newton's methods and the resulting linear system 
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can be solved very effectively by using the block-elimination method. A brief 
description of the solution procedure is given below. For details, see refer- 
ences 11 , 12 and 13. 

Let us consider the system given by (33) and (35). Using Newton's method, 
we can write the linearized difference equations for (33) as: 


h,. 


5f 




^i 


(47a) 

(47b) 


(s,)j6Vj + (S2)j4vj + (S3)j6f. + (S4)j6fj_i + (S5)jiuj + (s6)jSUj_, = (r^). 

(47c) 

Here 


= Vj-i/2 


-n-1 




n / 2 \ n 


(48a) 

(48b) 

(48c) 


/_ \ _1 , loO 

- HT+ r-fj 

vJ 

(=2>j 2^fj-l 

J 

pn 

pO 

(c \ ~ 1 wH-1 

^ 4^j ~ 2 ^3 


(49a) 

(49b) 

(49c) 

(49d) 


(S5)j = -P%n_^ 


-i 


(49e) 

(49f) 
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The boundary conditions given by (35) become 


STq = 0 , 6Uq “ 0 , 6Uj = 0 


( 50 ) 


The system (47), (50) has a block tridiagonal structure. This is not 
obvious and to clarify the solution we write the system in matrix-vector form. 
We first define the three-dimensional vectors 6. and r.- for each value of 
j by 








~ 0 

6* = 


0 £ j _< J 


0 


L J 











1 £ j 1 J-1 

= 



. 



0 


(51) 


and the 3x3 matrices A-, B*, C- by 

J J J 



1 

0 

0 


1 

-hj/2 

0 



> 

o 

II 

0 

1 

0 


(= 3)0 

(=5)j 


1 j 

£ J-1 


0 

-1 

-h^/2 


0 

-1 





1 

-hj/2 

0 


-1 

-hj/2 

0 



> 

C_1 

III 





(S4)j 


(= 2 >j 

1 1 J 

< J (52) 


0 

1 

0 


0 

0 

0 




0 

0 

0 








0 

0 

0 





0 < j 

< J~1 


0 

1 
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In terms of the above definitions, it can be shown that the system (47) and 
(50) can be written as 


where 







/A = 



Rj -1 Aj_i Cj^i 


— - 


“• — 

^0 


^0 

'^1 
r\j 1 



r\,c. 



• 


• 

• 

r. = 

'X.J 

* 

• 



• 


. 

^ J -1 


^ J -1 



P _ 


( 53 ) 


(54a) 


(54b) 


The solution of (53) is obtained by the procedure described in reference 13. 

Similarly, for the general case, if we assume (fj 

v"-’-’) and (f'?-’’’'. u';-’’’. vH-’’’) to be known for 
J J J J J J 

0 ^ j £ J, then the differenced equations for the general case can be solved by 

a similar procedure used for the two equations for the initial conditions. 
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IV. RESULTS 


The present method is applicable to both laminar and turbulent boundary 
layers. The flow is laminar at the leading edge, and it becomes turbulent at 
any specified x > 0. The initial conditions along the (n*t) plane are gen- 
erated by solving the system given by (30) and (35), and the initial conditions 
along the (n,x) plane are generated by solving the steady- flow equations, that 
is, (27) subject to (25). 

The marching procedure is along the t-di recti on. For a specified 
x-location, the governing equations are solved for each specified t-station. 
Since the linearized form of the equations are being solved, we iterate at 
each t-station until some convergence criterion is satisfied. For both 
laminar and turbulent flows we use the wall -shear parameter f^ as the con- 
vergence criterion. Calculations are stopped when 

|6f;| < 6^ (55) 

where the value of 6-| is prescribed. 


4.1 Laminar Flows 

There are several problems in which it is necessary to account for the 
fluctuations in the external flow. These fluctuations may change both in 
direction and in magnitude. A simpler case is one in which the external flow 
fluctuates only in magnitude and not in direction. This problem has been 
studied by Lighthill for a flat-plate flow. According to his analysis, 
for an external flow in the form 

Ue(t) = Uq( 1 + B cos o)t) (56) 

the reduced skin-friction coefficient c^/2 is given by two separate 
formulas depending on whether the reduced frequency cox/u^ is much smaller 

or much greater than one. 


"f 

T~ 



0.332 + 8(0.498 cos cot — 0.849 ux/u^ sin wt), wx/u^ « 1 
0.332 + B(o)X/Uq)^'^^ cos (wt + tt/4), cox/u^ » 1 


(57a) 

C57b) 
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In many problems it is often desirable to find the phase angle between the 
external flow and say the reduced skin-friction coefficient, which in terms of 
the transformed variables defined by (20) and (21), is f^. To determine this 
numerically for a fixed x = x^, let us consider the general case in which the 
external flow is given by 

Ug(x,t) = Uq(x)(1 + B cos wt) (58) 

and let 

f^(x,t) = g(x,t) (59) 

We use the following procedure to determine the phase angle between Ug(xQ»t) 
and g(xQ,t). We first compute lIq(x^) and gT(xQ) from 

to+P 



where P = 2 tt/o) and is the period of oscillation of the maintstream. Strictly 
Uo and g^ are functions of t^ as well as of x^, but for B 0.150, it 
was found that transient and other effects have virtually disappeared when 
t^ = P and this value of t^ was chosen to compute them. 

From (58 )» we can write 

where A = Uq(x^)B. 

Similarly, we can write 

g(xQ,t) — ^(Xq) = C cos [o)t + 4 i(Xq)] = C [cos wt cos <}>(Xq) — sin ojt sin 4 >(Xq)] 

(62) 

with (j)(xQ) denoting the phase angle between u^ and g at x = x^. If we 
take the product of (61) and (62) and integrate the resulting expression, we 
find cos 4>(Xg) to be given by 
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to+P 

f -g(Xo)]}dt 

cos +(x„) =.-^2 (63) 


Here 



[Ue(Xo,t) - UQ(xQ)]^dt 



to+P 

J [g(xQ,t) -?(Xo)]^dt 


(64a) 


(64b) 


Figure 4 shows the computed phase angle (between wall shear and external 

velocity) variation with reduced frequency wx/Uq for a laminar flat-plate 

flow. These calculations were made for an external flow given by (56) with 

B = 0.150, u = 17.5. A total of 41 t-stations, and 26 x-stations with 
0 

At = 0.20 and Ax = 0.1 were taken. Initially n was taken as 8 with 

• ''00 

An = 0.28. During the calculations n^ was allowed to grow. Our results 

(not shown in the figure) indicate that it is desirable to compute for two 

periods in order to determine the phase angle accurately. Figure 2 also 

shows the results obtained by Lighthil 1 with low frequency and high 

frequency approximations. These numerical calculations are in good agree- 

f 15l 

ment with those computed by Ackerberg- and Phillips^ they show that 

if B £0.150, Lighthill's low-frequency approximations to the phase angle 
are accurate for w £0.2u^/x and his high frequency approximation is 
accurate for u _> 2.6u^/x. 

If we write the wall shear as 

T " To (wt + (J))* 


then it follows that 


Zi = 

Tq ~ '1^ cos (o)t + (j)) 


20 




Figure 4. Variation of phase angle between wall shear and oscillating 
external velocity for a laminar flow over a flat plate. 

where is the wall shear for Blasius flow, then again t-j is strictly a 
function of x, B and t. However, with defined by (63) it appears that 
T-| in independent of B and the t-dependence has virtually disappeared when 
t > P in the calculations. It, thus, appears that nonlinear effects are 
negligible for B £0.150. In terms of our transformed variables . 


-1 _ 

Tq cos(iot + (f) 


( 66 ) 


and Figure 5 shows the computed values of for t = 2P together with 

Lighthill's low frequency and high frequency approximations. The agreement is 
good, the error at large tox/u^ being approximately equal to ^ in magnitude. 


4.2 Turbulent Flows 


The experimental data on unsteady turbulent boundary layers is very limited. 
The only data known to the author is that of Karl sson^^^^ , which consists of an 
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6.0 



figure 5. Variation of amplitude of the wall shear fluctuations for an oscillating laminar flow over a 
flat plate. 



oscillating free stream in a zero pressure gradient turbulent flow. The 

3 

measurements were made for a Reynolds number of R^*(=Uq6*/v) = 3.6 x 10 . 

H^re 


.../(' -y* 

In order to compare our numerical calculations with this data, we take the 
external flow to be the same as (56) and for convenience write it as 

'V/ 

Ue = Uq + u^B cos tot (67) 

(1 ) 

Here u^B = uj^ ‘ according to Karlsson's notation. 


We follow Karlsson's notation and denote the x-component of the velocity 
within the boundary layer by 


u(x,y,t) = u(x,y) + u^^^ cos cos tot — u^^^ sin ^ sin tot 
But since f = u/u^, we can write the above expression as 

Uo^' = ^Cx,y) + u^^^ cos (j) cos ut — u^^^ sin sin tot 
Using (67), we write (69) as 


(68) 


(69) 


, 0 ) 


u _ f _ u . u^ ^ cos (f) 4 . U'" sin 

— ^-yy - Tt; — - 7 ^ + |7(T1 ^ ^ SIH (OC 


B u„B 


,( 1 ) 
TTI 


(70) 


Karlsson gives lT/u^, u^^^ cos <|)/uf,^ ^ and u^^^ sin for differ- 

ent frequencies (to/Zir) and for different free-stream amplitudes B. To 
compute u^^^ cos <l>/up ^ from (70) we multiply both sides by cos tot and 
integrate with respect to t from 0 to 2co/'rr to get 


. 0 ) 


27r/to 


fl/ \ _ u^’' COS <|) _ a 1 f XI/ 4.\ 

A(n) = y V J f (n,t) COS 


totdt 


( 71 ) 


Similarly, to get u^^^ sin <i)/u^^ , we multiply both sides by sin cot and 
integrate. That gives. 
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/^x 2Tr/to 

B(n) 5 ( 1 ° - 7 1 f f (n,t) sin utdt 


( 72 ) 


To compute A(ri) an(d B(n) we take uniform steps in At and noting that 


A+ - 


t = nAt 


n 


n Ij 2j N 
n _ 


we approximate (71) and (72) by [recalling (t')j ~ ^ (nj»tj^)] 


(73) 


BN ^ 


n=0 


C74) 


N 

B. = (fO*^ sin o)t 

J RN ^ j " 

n =0 


(75) 


Figures 6 to 8 show the results for Karlsson's data. These calculations 
were started as laminar at x = 0 for the external flow given by (67) with 



Figure 6 , Comparison of calculated and experimental velocity prof;[les at 

R 5 * E 3.6 X 10^. Symbols denote u/u for values of B = 0.292, 
0.202 and 0.147. 
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^PRESENT METHOD 
ttNASH, ET AL 


Figure 7. Variation of phase angle between wall shear and oscillating external 
velocity for laminar and turbulent flows over a flat plate. The 
symbols denote Telionis' interpretation of Karlsson's data. 




Figure 8. Comparison of calculated and experimental values of the in-phase (A) 
and out-of-phase (B) components of an oscillating turbulent flow over 
a. flat plate, (a) Experimental data are for w/Ztt = 0.33 cycles/sec; 
ui^Vu^ = 29.2% (o), 20.2% (a), 14,7% (□). Calculations for for 
uP)/Uq = 14.7% 



y (cm) 

Figure 8. (b) Experimental data are for oj/27r = 1.0 cycles/sec: ui^ Vu = 35.2% 

(o), 19.5 %(a), 10,7% (□). Calculations are for uOJ/u^ = 35,2%. 
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(c) 


Figure 8. 



y (cm) 

(c) Experimental data are for u/Z-rr = 4.0 cycles/s 
(o), 13.6% (a), 6.2% (□). Calculations are for 




= 26.4% 
/Uq“= 26?4%. 


Uq = 5.33 m/sec. The Ax~spacing for turbulent flows was one foot. The turbu- 
lent flow calculations were started at x = 3.05 cm, and at z ^ 3.7 rn the 
experimental velocity profile was matched (see figure 6). We note that the 
experimental velocity profile is essentially unchanged for various magnitudes 
(small) of external flow. 


Figure 7 shows the variation of phase angle (between wall shear and edge 
velocity) with reduced frequency wx/u^ for laminar and turbulent flows. The 
results indicated by the present method were obtained by the procedure discussed 
above for a)/2TT = 0.33, 1, 4 cycles/sec with B = 0.147, 0.352 and 0.264, 
respectively. The other results shown in this figure are obtained from a fig- 
ure given by McCroskey^^^^ . In that figure, McCroskey gives a comparison 
between numerical studies by McCroskey and Philippe, Tel ion is and Tsahalis, 
and Nash et al . , and also some experimental data due to Karlsson interpreted 
by Telionis. The phase angle distribution for laminar flows is computed by 
the present method and is included for comparison purposes. 


Figure 8 shows the comparison between calculated and experimental values 

of the in-phase and out-ofrphase components, that is, cos 4>/u^^^ and 

( 1 ) ( 1 ) ” 

u' 'sin (()/u^ , respectively, for the three different frequencies and ampli- 

tudes mentioned earlier. It is interesting to note that of these three cases. 
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only the ^ — distribution corresponding to w/2Tr =» 1 has posi- 

tive and negative values, just like the computed onesf 

Table 1 presents the ratio T-j/T^ for laminar and turbulent flows. It 
can be shown that '^i/'To calculated from 


1L= ~ '' 

^0 6 cos (wt + (})) 


(76) 


Here cf^ denotes the flat-plate skin friction; it can be calculated from the 
following formula given by Cebeci and Smith ^ ' 




0 (2.44 In Rg* + 4.43)' 


(77) 


The results shown in Table 1 indicate that in turbulent flows the ratio tt/t 

I 0 

is much smaller than its value in laminar flows. We also observe that for 
turbulent flows the variation of with reduced frequency is relatively 

small compared to its variation for laminar flows; this is consistent with the 
results reported by McCroskey and Philippe^^^. 

Table 1. Variation of t^/t^ with Reduced Frequency, w 


'X, 

0) 

Laminar 

Turbulent 

1.36 

3.3 

2.28 

1.48 

3.50 

2.39 

4.50 

6.50 

2.75 
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V. CONCLUDING REMARKS 


A numerical method for computing two-dimensional incompressible unsteady 
laminar and turbulent boundary layers with fluctuations in external velocity 
is presented. The Reynolds shear stress term is modeled by using the alge- 
braic eddy-viscosity formulas developed by the author. The numerical accuracy 
of the method is first checked for laminar flows and the results are compared 
by Lighthill's analytical solutions as well as by those obtained by other numer- 
ical methods. Based on these comparisons it is found that the method is quite 
accurate. 

Similar calculations are also made for turbulent flows and the results 
are compared with available experimental data and with other numerical solu- 
tions. In general the present results agree reasonably well with experimental 
data and with those obtained by other numerical methods, McGroskey and 
Philippe^^^ and Telionis^^^ which also used the present author's eddy- 
viscosity formulation, Cebeci^^^. The present eddy-viscosity formulation 
isrslightly different than that reported in Cebeci^^^ in that it allows the 
unsteady effects through p^ in equation (18). Furthermore, it also accounts 

for the low Reynolds number effect which was not included in Telionis' calcu- 
( 2 ) 

lations' * . Another difference between the present results and with those 
obtained by McCroskey and Philippe and Telionis and Tashalis may be due to 
the procedure used in computing phase angles and the procedure used to gener- 
ate the initial conditions required to make the turbulent flow calculations. 

Although our turbulent flow calculations agree reasonably well with 
Karlsson's experimental data and with those methods that use algebraic eddy- 
viscosity formulas, they do not agree with those computed by Nash et al . 

This is somewhat surprising because that method is a modified version of 
Bradshaw's method, and that the predictions of this method for steady flows, 
see ref. (5), and unsteady flows, see ref. (10), agreed well with the pre- 
dictions of the author's method. 
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